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abstract 

The Advanced Sensor Development Labor = ( ASDUa, ,he Stenni^ Space Cemer tUve.ops, 
maintains and calibrates remote sensi g trade-offs analysis, and establish system 

Admin, strat ion (N ASA). To pe. * s ^f^^aWicaUimXL of sensor systems. This 
parameters, ASDL has developed a : soft' l nfraRpd Fneineerine " - ATTIRE, simulates the 
package called " Analytical Tools for of the sensor to be 

various components of a sensor system. performance parameters are then integrated to 

analyzed independently for its performance P (SNR) Noise Equivalent Radiance 

obtain system level informal, on such as S tadro „ desc 3bes the uses of the 

tNERl Noise Equivalent Temperature Difference (Nb 1 D) etc. i ms papci u 
package and the physics that were used to denve the performance parameters. 

The package runs from a menu-tlnven interface 

Se various'aspects of^' interactive fe’atums allow for 

data plotting and visualization. 

The package can simulate a •^-^“ZShelCon 

changing the appropriate parameters ^ versus any system parameter. A parameter (such as 

package can proyide plots of varied and EffccTon mother parameter (e.g., NETD) 

and then projects them onto a common wavelength array for all the paramete . 
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INTRODUCTION 


a^cls of sdcncc ^deng?^ring, ^rom'thc^un^rst^^^^^heh^s^f ? “ ° f ; *»«■ varied 
to the computation of signal-to-ratios of a spncnr th* • g f ^ c aws °f electromagnetic radiation 
suited for compute? s Sation In Ms **“? ""W* and desi ^ n is idea “y 

called ATTIRE runs on a PC ™OS 15 deSCribed ' The P ac W 
understanding the behavior and interaction of various sub-svstems rhar systems as as 

shows the flow-chart for the energy flow that is modeled bv ATTIRF Thf» nS ' UtC - 3 sensor ' Tig. 1 
its various parameters as shown. moaeied by ATTIRE. The sensor is represented by 


MENU DRIVEN SYSTEM 

-msrv 

b °,& f r a nux source ™° 

parameters that make up the sens!>rTux source* On dSmr^.h,. V™" ' S USed t0 rcprcsem ‘ he 
window which displays the characteristics of the parameters th£ il ? . one ent 5 rs int0 another 

flux energy is then propagated through ?he nex^bsvlif " P ** s u ource ^- The source 
attenuates the radiation coming from the source and also aSTn the atmos P] ,er e- The atmosphere 
to effects such as scanering. Th.s is due 

of the source onto a detector. The detector converts *Sado??nm elecoSf ene £ *£? 
system and the electro-optics ie the dpipnnr ct.h cc into electrical energy. The optics sub- 

The electrical energy can then be amplified and modified inrir^h 6 "h X -r!'h° IC ? nS the Signal flow - 
This is represented by the electronics icon The ~ bandwidth and noise characteristics. 

system, i.e., its dwell time, IFOV etc. The user can click ontoeSJ!?of?h Spat “ charactenstlcs of the 
sub-systems. me user cm click ° n to either of these icons to enter into these 

Apart from these six icons the user encounters the main menu Th» m „- 

window. The window can be resized or relS usme a m^i m me " U ,S ™ de U P of a 

edges of the window for resizing. The window ca/alsoTml?, ? pressed at any of the 

^a S re?^ The ouUine of ^ wmdo w appear^ SS 

channel °" a channel-by. 

The listing in the window contains information on the current sensor being analv^H Th* r ,• • 

done on a channe -by-channel basis Each channel ic k • u # analyzed. The listing is 

the main window shows system level ^ by its bandwidth. For each channel, 

parammen' rNrN E rand ,h N a ^ nC ' den ' ™ 

when the number ofrtwircl^ This happens 

‘CHAN* in the horizontal menu UmbefS P r °g ramm able by using the item called 

The horizontal menu performs functions to modify or analyze the sensor parameters. The first menu 
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item ‘CHAN’ performs three major functions. It allows the user to modify the number of channels in 
the sensor being analyzed. The maximum number of channels that currently can be processed is 50. 
The second function of this menu item is to modify the bandwidths of the channels. 

The process of modification of any parameter in ATTIRE follows a standard procedure. The user is 
first prompted with a menu asking if the modification needs to be done for only one channel, a few 
channels, or for all channels. For the case of one or a few channels, the user is prompted for the 
number of the channels that need to be modified. A tick mark appears next to the number of the 
channel. In the case of more than one channel being prompted for, the user can continue to click on 
the channel numbers. To delete a number that has been tick marked, the user just needs to click on it 
again. One continues this process until the required channels have been selected. The user then clicks 
outside the slection windows to accept the selected channels. Also, all the operations in ATTIRE can 
be done without a mouse control. In the case of channel selections, the up and down arrow keys are 
used to move within the selected window. If the number of channels is greater than the size of the 
selected window, the “Page Up” and “Page Down” keys allow for scrolling to other areas of the index. 
The tick marks are created by typing the “Enter” key. The end of the selection process is signified by 
the “Control - Enter” keys. After having selected the number of channels that need to be modified, the 
appropriate parameter is edited. 

If only one channel is selected, the existing value of the parameter appears in a window. This value is 
appropriately edited. If more than one channel is selected, the modification can be done in two ways. 
The first method just allows for placing the same parameter value in all channels. This procedure is the 
same as for one-channel selection. The other manner in which the selected parameter can be modified 
is by placing a linearly increasing series of values in the selected channels. In this case the user is 
prompted for the initial and final values of the parameter starting from the first to the last selected 
channel. 

The source emissivity, atmospheric transmission, and optical transmission are initially assumed 
constant over the specified spectral bandwidths. There is an option for making these parameters 
spectrally variable. If more than one of the above-mentioned three parameters are spectrally variable, 
then it is possible that the upper and lower wavelength values and the resolution of the wavelength 
array may not be the same for all three parameters. The software package performs an interpolation of 
the data to smooth out the curves and then projects them onto a common wavelength array. 

The spectrally distributed incident flux can be analyzed as it propagates through the subsystems. Fig. 1 
shows the signal flow and the various parameters that affect the performance of the system. Each of 
the icons on the main menu represent one of these sub-systems. 

SUBSYSTEMS 

The software divided the system into the following sub-systems: 'SOURCE FLUX', 
ATMOSPHERICS', OPTICS’, 'ELECTRO_OPTICS', ’ELECTRONICS’ and PIXEL 
PARAMETERS’, 

The system can be simulated in a multi- spectral manner, i.e as many as 50 channels can be processed 
simultaneously for performance analysis. For each channel, the parameters corresponding to these sub- 
systems can be entered. The parameters can be entered one-by-one for each channel or in a linearly 
varying manner to analyze the effect of varying one parameter on the system performance. Also, all 
the channels can be made to assume the same value. 

Since each channel is modeled independently of the other channels, the package can be used to 
simulate one 50-channel sensor or 50 unique sensors or any combination in between. 
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FLUX SOURCE 


For analysis of the system a target has to be simulated as a source for the flux signal. This could be 
either a black body for thermal sources, the solar irradiance for visible/near IR sources or any 
customized source depending upon the application. Upon entering the menu the user can choose either 
of these options. 

Blackbodv Source 

The spectral distribution of blackbody radiant exitance in energy units at a given temperature is 
obtained from Planck's formula 

F n (T) = f M eX a, T) <R n (X)dX 
Jax„ ’ 

where 


AX„ = spectral bandwidth of channel 'n'. 

SK n (X) = Relative Responsivity of the sensor for channel 'n'. 


M e T) = ^ W cm -2 pm -1 

X 5 ( e xt - 1 ) 

where 

c, = 27thc 2 = 3.7483 x 10 4 W cm -2 pm 4 

c 2 = — = 1.4388 x 10 4 pm K 

h = Planck's constant, 
k = Boltzmann constant. 

X = wavelength (pm). 

T = temperature of source (K). 


( 1 ) 
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Assuming the source to be Lambertian , the spectral radiance is given by dividing the exitance by n 
steradians. If the source of the flux is not a blackbody, the spectral radiance is obtained by multiplying 
by the emissivity of the source as follows 


4jl<VT) = ^ 


eft.) 


71 


'\ 


, C 2 

V X 5 ( e XT - 1 ) 


_2 -l 

W cm Jim sr 


where e(X) = emissivity of the source. 


( 2 ) 


For a field-of-view (FOV) ft s , the radiant exitance is given as follows: 

M e>x a,T) = (Q s ) L e X aj) Wcm" 2 pm" 1 

where Q s = Field of View. (3) 


Solar Source 

Solar radiation can be approximated by a 6,000 K blackbody curve. The solar spectrum is modified 
by atmospheric transmission as it passes through an air mass en route to the Earth's surface. For an 
airborne remote sensing system, this energy is further affected by the albedo of the surface type and 
the atmospheric path back to the sensor system. The spectral distribution of the solar irradiance was 
compiled from various sources 12 to arrive at a curve for wavelengths ranging from 0.295 to 2.541 
pm. The resolution of this curve is 1 nm. This is the solar curve incorporated in the ATTIRE 
program, and used to calculate the available solar radiation for each reflecting channel. Fig. 3 shows 
the solar curve. 

Modeling of the Source Flux in ATTIRE 

ATTIRE provides an icon in the main menu called 'FLUX'. Fig. 4 shows the flux menu as it appears 
on the screen. On entering this menu from the main menu, another window appears along with a 
series of icons. These icons allow the user to change the type of flux ( Blackbody or Solar ), the 
temperature of the radiating source, the emissivity or reflectivity of the source and the Field-of-View 
(FOV) of the source. The window contains a horizontal menu whose items allow the user to perform 
the same operations as the main menu. 

The flux type can be changed by clicking on the icon or typing 't' key. On entering the icon, the user 
is prompted for identifying the type of the flux. If the type is solar, care should be taken that the 
bandwidths for that channel are within the solar range. The Solar menu contains items that identify the 
solar flux; i.e., channel bandwidth and reflectivity of the target. The user enters the appropriate 
bandwidth values and the program picks up the corresponding irradiance values from a file called 
SOLAR.BIN. It then multiplies each irradiance value by the reflectivity to obtain the effective spectral 
flux. Selecting the Other Sources menu, the user can enter the source as a specific lamp or any 
emitting source whose distribution is known. 

To obtain spectral radiant exitance for a source radiating uniformly in all directions, the FOV is made n 
steradians. To obtain spectral radiant exitance at any other solid angle, the appropriate value is entered. 
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In order to compute radiance, the FOV should be set equal to one steradian or 57 T The FOV is 
modified by clicking on the icon or by typing V key. 


If the user has selected a thermal source as the flux, the program uses the temperature assigned for that 
channel to compute the flux. The temperature icon is used (or typing ‘t’) for that purpose. 

The medium between the source and the detector is made up of the atmospheric path and the optical 
path. The flux passes through the atmosphere before being collected by the optics of the sensor It 
then passes through a senes of reflective and refractive optical elements before it is incident on ‘the 


Depending on the type of flux, the radiation is attenuated by the emissivity for the thermal and the 
albedo for the solar case. These are modified by clicking on the appropriate menu (or by typing ‘e’). 
This value can be spectrally variable or constant. In the case of spectrally variable values, the user is 
prompted fora file name. This file name must contain the two arrays - wavelength and emissivity or 
reflectivity. This is an ASCII file which has these two parameters in each line of the file. The first two 
fines of the file are ignored and can be used for comments. The data is picked up from the third line 
The first number is the wavelength and the second the parameter. The file continues to read until the 
line it read does not contain two numbers or end-of-file is reached. When a spectrally variable 
parameter is chosen, the parameters entered in channel bandwidth are not applicable as a wavelength 
array supplied with the file determines that bandwidth. 6 


ATMOSPHERE 


The atmosphere attenuates the amount of flux that enters the sensor. Only a fraction of the flux is 

incident on the optical aperture. This atmospheric transmission coefficient is denoted by x (A.) It 

is assumed to be a function of the wavelength. a 


Also, the atmosphere contributes to the radiance entering the sensor system. This atmospheric path 
radiance can be divided into two categories, depending on the direction of the flux. The upwelling 

radiance, LA' (A.) , is the flux from atmospheric particles that is directly incident on the sensor. The 


downwelling radiance,LA2 (A.) , is the flux that is incident on the target. This radiance then is reflected 
rrom the target back to the sensor. 


Thus the effective flux incident on the sensor is given by: 

L*. (A.,T) = x a (A.) [ e (A.) lJ? a,T) + l£‘(A.) + (l-e(A.)) l£ 2 (A.) ] 
A H IRE Modeling of the Atmosphere 


(4) 


TTie atmosphere is selected by clicking on the icon for the ‘ATMOSPHERE’ in the main menu. Fig. 5 
shows the menu as it appears on the screen. This menu contains three items for each channei - 
transmittance upwelling and downwelling path radiance. The transmittance is modified in a manner 
similar to that for the emissivity in the flux sub-menu. 


OPTICS 


The aim of the optics subsystem is to collect, focus, and disperse the radiant flux from the source. The 
energy is then focused on the individual detectors. The two parameters that are critical to the collection 


506 


of this energy are the area of the collecting optics and the focal length. 

The radiance, limited by the solid angle subtended by the ground pixel, is incident on the entrance 
aperture and determines the irradiance available to the optics. The area of the entrance aperture minus 
any obscuration is the other limiting factor. 

The solid angle is given by 



where A pix = area of the ground pixel 

a = distance between the source and the entrance aperture (altitude) 

This solid angle must equal the solid angle subtended by the detector on the entrance aperture 
(invariance theorem). 


A A, 



where A d = area of detector 

efl = effective focal length ^ 

Depending on the design of the optical system, the entire area of the primary mirror may not be 
collecting the incident energy. This could be due to an obscuration of the mirror. 

Aeff = Area of primary mirror - Area of obscuration 

The series of reflective and refractive elements that constitute the optics subsystem are modeled by an 
optical transmission coefficient x Q (X) , the effective area of the collecting optics (A e ff), and the 
effective focal length (efl). 


In the case of several systems , the energy is focused on the field stop so that the beam can be directed 
to more than one detector. For such systems, the area of the field stop should be substituted tor A d . 

ATTIRE Modeling of the Optics 

The main menu of ATTIRE contains an icon for the Optics. On clicking this icon, the user encounters 
a window that contains the optics parameters. Fig. 6 shows the menu as it appears on the screen. 
These are the diameter of the entrance aperture, the diameter of the obscuration, the effective local 
length, and the optical transmittance. 

The diameters are modified by entering into the appropriate menu items. The effective focal length is 
modified through its menu item. This value is also dependent on the area of the detector and the IFOV 
of the system. For all the flux to be incident on the detector, the following relationship must hold. 
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2 


where a = EFOV 


(7) 


Hence, by changing the focal length, the program automatically changes the detector area according to 
the above equation. 

DETECTORS 


For a background noise limite d infrared ph otoconductor (BLIP), the theoretical spectral D* is given by 
D*(A.) = — : — cm Hz 1/2 W _1 


2hc sin 0 1/2 ^bb ^g) 


where 

T) = Quantum efficiency of the detector 

0j/2 = Cold shield half angle (radians) 

BB 

Mp = Photon background noise flux 
= JpkM®?a,T B m 
T b = Background temperature (°K) 

( 8 ) 


The peak D* of the detector is obtained by substituting for peak wavelength in the above equation. The 
D* curve for the detector is now created by linearly extrapolating D*( ) up to the peak wavelength. 

D* (X) = D* (JL k ) cm Hz 1/2 W _1 

^ P k (9) 

D** is D* normalized to cold shield half angle and is given by: 

D** = D* sin 0 1/2 (1Q) 

A1T1RE modeling of the Electro-Optics 

The main menu of ATTIRE contains an icon for the detectors. On clicking this icon the user enters a 
screen that contains a window for the parameters that make up the detector sub-system. Fig. 7 shows 
the menu as it appears on the screen. There is an accompanying icon for each parameter which allows 
the user to change their values. The user can use the theoretical model to simulate the peak D*. Also, 
a user-defined peak D* can be input. If a peak D* is input, the program automatically adjusts the other 
parameters, e.g., Quantum Efficiency, to provide the input D* peak value. The D* spectral curve is 

created by drawing a straight line from X=0 (where D*=0) to X pk (where D* pk is given). The 
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detector area does not change independently. As discussed, it is dependent on the effective focal 
length of the optics and the IFOV of the system ( Eqn. 7). 

PIXEL PARAMETERS 


This sub-system models the the pixel geometry as defined by the sensor motion and system spatial 
IFOV. In the case of airborne scanners, the sensor is mounted on an aircraft, whose velocity and 
altitude of flight need to be controlled for contiguous scanning. 

The pixel size is a function of the altitude (a) and IFOV, and can be computed by 


, a 

Pixel Size = 2 a tan[ y 


(ID 


where a = angular IFOV 

a = altitude of aircraft in meters 

Pixel size can also be determined using the small angle approximation by: 

Pixel size = (IFOV) (a) 

For contiguous scanning, the time required for one revolution of the scan mirror is equal to the time 
required for the aircraft to travel forward a distance equal to one IFOV on the ground. Using the small 
angle approximation for pixel size 


n 


where 


aa 

v 

n 

v 


( 12 ) 


= scan speed in rps 
= aircraft velocity in meters/sec 


ATTIRE Modeling of the Spatial Parameters 

The spatial parameters are modeled by entering the menu entitled 'PIXEL PARAMETERS . Fig. 8 
shows the menu as it appears on the screen. Since these parameters are interrelated, changing one may 
also affect the other. This menu contains the following items: 

SCAN-SPEED-HT, SCAN-SPEED-VEL, DWELL-TIME, PIXEL-SIZE, and IFOV. 

The scan speed is a function of both altitude and velocity (Eqn. 12). 

The dwell time is a function of the scan speed, which in turn is a function of the altitude and velocity. 
In ATTIRE, changing the dwell time affects the scan speed and the velocity only. 


l d - 


a 

2nn 


a 


- . . seconds (131 

2n (v/a) 

The Pixel menu contains the IFOV, pixel size, altitude, scan speed, velocity and dwell time. Several 
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of these parameters are inter-dependent. Therefore changing one might affect some of the others as 
well. 


ELECTRONICS 


The signal conditioning and preamplifier electronics are modeled for their bandwidth and noise 
properties. The aim of the sensor design is to be detector-noise limited. In practice, some noise is 
added by the electronics. Hence the total noise of the system is represented as follows: 


Noise total — Noisedetector Noise e lectronics 


The electronics noise is input as a noise factor. This is defined as the ratio of the total system noise to 
the noise of the system prior to the electronics. This makes it simple to define the noise as a multiple 
of the rest of the noise. Also, it indicates the contribution of the electronics to the total system noise. 
The detector noise is accounted for in the model where the D* is calculated assuming certain 
background noise. 

The noise of the electronics is defined as follows 


(SNR) input 

n f = 

(SNR) output (14) 

The electronics also determines the sampling interval required to sample the analog output of the 

amplifiers. In the case of aircraft scanners, as stated earlier, the scan speed ‘n’, the IFOV ‘a’, the 
plane velocity ‘v’ and the altitude ‘a’ are related as follows: 


a 

td - 2nn 


a 


2 k (v/a) 


seconds 


(13) 


The sampling period is defined by the dwell time of the sensor i.e the time the system spends at each 
pixel. In the case of a scanner this is given by 


The dwell time determines the smallest spatial sampling interval. It is the smallest time interval 
between which any occurring change can be detected, i.e., any change faster than td will not be 
detected. It is the time interval between two successive pixels. The spatial sampling frequency is the 
inverse of the dwell time. This means, based on the Nyquist criterion, the largest spatial frequency 
component present in the signal is one half of the spatial sampling frequency. Thus in order to collect 
this signal, the bandwidth of the signal conditioning electronics must, at most, be half the spatial 
sampling frequency. 
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Dwell Time Spatia ] Sampling frequency 
Bandwidth = j Spatial Sampling frequency. 


(15) 


ATTIRE modeling of the Electronics 

The user enters into the Electronics sub-system by clicking on the icon in the Main Menu. The noise 
factor and the Electronics Bandwidth can be modified independently for each channel. Fig. 9 shows 
the menu as it appears on the screen. 

COMPUTING THE BAND-PASS FLUX 


The spectral radiance incident on the detector (assuming no path radiance) is given by. 

, x c x z(X) xJX) T 0 (W -2 -1 

L e X (X,T) = — w cm sr 


n r (exp (C 2 AT) -1) 


The solid angle subtended by the source at the entrance aperture is 


— = Q c steradians 
2 s 

a 

Apjx = Area of the pixel 

= (aa) 2 

The total energy incident on the detector is given by 

M (e .n) (ti, T) = j ^ 2 M e> x Oi,T) dX = 2 Q s Ot,T) d>. W cm 


(16) 


(17) 


( 18 ) 


where X x and X 2 are the lower and upper wavelengths of the channel. 

For cases where the bandwidth of the channel is defined by a relative spectral responsivity curve, the 
total energy is given by: 


Jch# V n J X 5 (exp(C 2 AT)-l) 
where R(X) is the relative responsivity curve. 


(19) 
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Power Incident on the Detector 


The power incident on the detector is obtained by integrating the radiance exiting the imaging lens over 
the channel bandpass and multiplying the integral by the throughput. This is given by 


p d = ( J e (X) x a (X) t g (X) L x (X,T) dX ) (y d ) 
where y d , detector throughput is given by 


— A d ft d 


where 


A d = Area of the detector 

Q d = Solid angle subtended by the detector at lens 
_ Area of lens 

(efl 2 ) (20) 

Rewriting using the above relationships, we obtain 

p d = f f e a) x a (X) t 0 (X) l x a,T) cal (A d ) 

Uax ) (21) 


Noise Equivalent Detector Power 

The noise equivalent power (NEP) of the detector is defined as the power incident on the detector such 
that the signal-to-noise ratio is unity. The NEP can also be obtained if the broad-band D* is available. 

In this study, the D* at peak wavelength is modeled. The broad-band D* can be obtained from the 
peak spectral D* as follows: 


Dbb (T) = 


f D* (X) L e X a,T) dX 

•'O 

f L e>x a,T)dX 

•Jo 


The NEP is then defined as 


(22) 


NEP(T) = 


V A d Af 

D^b (t) 


where Af = electronic bandwidth 


( 23 ) 
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Signal-to-Noise Ratio 

The signal-to-noise ratio (SNR) of the detector can be obtained as follows: 


SNR d = 


Pd 

NEP 


(24) 


Noise-Eoui valent Radiance 


The noise equivalent radiance (NER) denotes the required bandpass radiance exiting from the target 
such that the SNR is equal to unity of the system. 


NER = 


Bandpass radiance exiting the target 
SNR 


!*( XJ) dX 

'AX 

SNR 


( 25 ) 


Noise-Equivalent Temper ature Difference 


The radiometric performance of a thermal sensor is determined by its NEAT^ This is the t ^ 
sensor to theoretically discriminate between two temperature values. It is defined as the tempe at 
difference required to make the signal-to-noise ratio of the sensor unity. 


^ AT 3T tiT 

NEAT- snr" _^V s j " ^ Vs 

where V s and V n are the signal and noise voltages respectively. 

v,. = (solid angle) ^(>.)L^(X, T) dX. 

Jx, 

v 

^X) = Spectral Resposivity = r D*(A.) 

V A 

Expanding the expression and substituting the appropriate parameters 


NEAT = ^ pk - 

d a pk ) 



efl 2> | 1 
A eff J (Integral) 


( 26 ) 


( 27 ) 


( 28 ) 
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where Integral - 



x a (X)x 0 (X) e(X)X 


d L® B a,T) 

dT 


dx 


(29) 


GRAPHICS 

The program provides extensive graphics capabilities for viewing data parameters The graphing 
routines are divided into two categories - PLOT and GRAPH. The GRAPH feature is used 
exclusively for spectral graphs, i.e., the x-axis is always the wavelength. In the main menu, the 
GRAPH feature allows for plotting of either the Energy or the Photon Flux versus wavelength. The 
flux can be plotted for either of the existing channels or for any user-defined set of parameters. The 
program prompts the user to select one of these two options. If the user desires a spectral plot of a 
channel, only the channel number is required. In the case of a user-defined graph, the user is 
P h°cnA/ d f ° r thC tyPC fluX ( thermal / solar )> die temperature (if thermal), the emissivity/albedo, and 


The graph feature in the FLUX, OPTICS, and ATMOSPHERICS sub-systems are used to graph the 
spectral emissivity, the optical and atmospheric transmittance respectively. These can be variable or 
constant depending on the channel. Each of these three sub-systems also allow for an option in their 
GRAPH Menu to graph all these three spectrally variable parameters on the same graph. 

The plot item in the main menu is used to plot non-spectral parameters against each other. One such 
example is to compute the SNR of the sensor as it varies with the temperature of the source. The user 
is prompted for the bandwidth of the sensor. All other parameters that enter into computing the SNR 
are used from a thermal channel. The user is also prompted for that thermal channel number Next 
the range of temperatures over which the data needs to be plotted is asked for. 

The other significant feature of the PLOT item is the User-Defined feature. This allows the user to 
select either of the sensor parameters as it x and y axes. Then for all the existing channels of the 
sensors the plot is created. 

The graphs appear in a resizeable windows form. One window displays the graph and the other 
window the contents of the array that are being graphed. So the user can not only view the graph, but 
also view the actual values being plotted. The graphs can be resized by pressing the mouse or the 
edges and sliding it to the required window size. The top right button of the graphics window is used 
to make the graphics window fill the entire screen. This can also be done by typing the T key (for full 
screen). On clicking the top right button of the array values window, the user can quit the graphics 
environment. This can also be done by typing the ‘q* key. The array values window can be scrolled 
up or down to view different parts of the array. 

The main menu additionally contains an item for Graph menu to plot the energy or photon flux of all 
the channels. It also allows the user to plot all the spectrally variable parameters - which might include 
emissivity and atmospheric & optical transmission for a channel - on the same screen. This item 
a lows for obtaining graphs of spectrally varying parameters. Other graphs can be obtained from the 
plot menu. 

An additional option exists for the user to be able to obtain a hardcopy of these graphs on an HP 
Laserjet II printer. Fig. 10 shows one such graph. These figures are the screen outputs of ATTIRE. 
This is the result of plotting the energy flux radiated at various temperatures within two prominent 
thermal bands - 3-5 and 8-12 microns. At low temperatures, the 8-12 channel has more energy. 


514 


However as the temperature of the target increases, the 3-5 channel has more energy. Thus to monitor 
high temperatures ( specifically at temperatures greater than 582K ), e.g., volcanoes and forest fires, it 
is better to use 3-5 micron channel and for low temperatures, e.g., normal earth temperatures, at 300 K 
the 8-12 micron channel gives a better performance. 

PLOT Y vs X 

The package can provide plots of system performance versus any system parameter. A parameter X 
e.g the entrance aperture of the optics, could be varied and its effect on another parameter Y e.g 
NETD, can be plotted. Another parameter Z e.g the obscuration, could be varied for each plot and 
several plots obtained on the same graph. The menu for such " Y vs X plots for different values of Z 
" contains various options in each sub-system. Fig. 1 1 shows the result of plotting the variation in 
NETD with temperature for three different thermal channels. 


SPECTRALLY VARIABLE PARAMETERS 

The emissivity, atmospheric transmission and the optical transmission default as constants over the 
specified spectral bandwidths. These could also be made variable. Indeed, in typical situations these 
are spectrally variable parameters. Before entering either of these parameters from their respective 
menus, the program asks the user if they are variable or a constant. If they are a constant, the number 
as a fraction of 1 is entered in the box. If it is a variable, the box shows the letters VAR. The software 
then picks up the spectrally variable curves from the disk and defaults the bandwidths for the channels 
to the corresponding bandwidth and resolution. In the case of atmospheric transmission, for each 
channel the curve should be stored in ASCII format as al.dat, a2.dat, ... a6.dat. Similarly for optical 
transmission the first letter should be 'o' ( ol.dat .. ) and for the emissivity it should be 'e' ( el.dat .. 
). The first two lines of the data files are ignored by the reading routine. These shall be for the user to 
describe the curve for their purposes. The data shall follow from the third line onwards. Each line 
contains the wavelength ( in microns ) and the corresponding parameter ( in % ). Thus the reading 
routine picks up the values for the wavelength and the parameter, one at a time from each line until it 
reaches the end of the file. 

Interpolation and Projection 

If more than one of the above mentioned three parameters are spectrally variable, then it is possible that 
the upper and lower wavelength values as well as the resolution of the wavelength array may not be the 
same for any two out of the three arrays. For example, the emissivity data may be available from 3.0 
to 5.0 microns in intervals of . 1 micron whereas the optical transmission data may be available from 
2.2 to 4.85 microns in intervals of .05 microns. ATTIRE contains a methods for resolving such 
ambiguities. 

The package first checks for these discrepancies before deciding on the least common denominator for 
the resolution and selecting the lower and upper limits for wavelength such that all the existing points 
available can be plotted as they were obtained. However doing this results in several gaps in the 
arrays, e.g in the above example, a resolution of .05 microns and a bandwidth of 2.2 to 5.0 microns is 
chosen. This leaves blanks in the emissivity array at 2.2 - 3.0 microns as well as at 3.05, 3.15 ... 
microns. 

The package performs an interpolation on the data to smooth out the curves and then projects them 
onto a common wavelength array for all the parameters. The emissivity from 2.2 to 3 microns is not 
assumed to be zero throughout but starting from the value at 3.0 microns gradually tends to zero at 2.2 
microns. Similarly the emissivity at 3.05 microns is the average of the emissivity at 3.0 and 3.1 
microns. 
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SUMMARY 


The design of a visible through thermal IR sensor system requires a detailed analysis of how the input 
signal propagates through the system. The major components in developing a model for the sensor 
system are the source, atmosphere, optics, detector, spatial parameters, and preamplifier 
electronics. The final goal of the analysis is to determine the NER for the various spectral channels 
of a sensor system. 

In this paper, a simulation package "ATTIRE", for analyzing sensor systems was introduced. The 
package runs in a PC-DOS environment and consists of one executable program and several 
supporting files. The entire package fits on one high-density floppy disk. 

ATTIRE is a useful tool for performing design trade offs as it inter relates several aspects of the sensor 
system to yield performance parameters. It is also useful as a tool for the understanding of the 
concepts of radiometry. 
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